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Abstract 

In the framework of the littlest Higgs model with T-parity (LHT), we study the 
production processes of T-even (T + ) and T-odd (T_) partners of the top quark at the 
Large Hadron Collider (LHC). We show that the signal events can be distinguished 
from the standard-model backgrounds, and that information about mass and mixing 
parameters of the top partners can be measured with relatively good accuracies. With 
the measurements of these parameters, we show that a non-trivial test of the LHT can 
be performed. We also discuss a possibility to reconstruct the thermal relic density of 
the lightest T-odd particle Ah using the LHC results, and show that the scenario where 
Ah becomes dark matter may be checked. 



1 Introduction 



The hierarchy problem in the standard model (SM) is expected to give a clue to explore 
physics beyond the SM. This problem is essentially related to quadratically divergent cor- 
rections to the Higgs boson mass, and it strongly suggests the existence of new physics at 
the TeV scale. At the new physics scale, the problem is expected to be resolved due to the 
appearance of a new symmetry which controls the Higgs boson mass. With this philosophy, 
a lot of scenarios have been proposed so far. The most famous example is the supersym- 
mety (SUSY), by which quadratically divergent corrections to the Higgs boson mass are 
completely cancelled. Another example is the Gauge-Higgs unification, by which the gauge 
invariance in higher dimensional space-time protects the Higgs potential from any ultraviolet 
(UV) divergent corrections. 

In this article, we consider the third possibility, so-called the little Higgs (LH) scenario [TJ, 
in which the Higgs boson mass is controlled by a global symmetry. In this scenario, the Higgs 
boson is regarded as a pseudo Nambu-Goldstone boson arising from the spontaneous breaking 
of a symmetry. Due to the symmetry imposed, new particles such as heavy gauge bosons 
and top-partners are necessarily introduced, and main quadratically divergent corrections to 
the Higgs boson mass vanish at one-loop level due to contributions of these particles. Unlike 
the SUSY scenario, the cancellation of the quadratic divergence is achieved only at one-loop 
level, thus the LH model needs a UV completion at some higher scale. However, due to the 
cancellation at one-loop level, the fine-tuning of the Higgs boson mass is avoided even if the 
cutoff scale of the LH model is around 10 TeV. As a result, the LH model solves the little 
hierarchy problem [2]. 

Unfortunately, the original LH model is severely constrained by electroweak precision 
measurements due to direct couplings among a new heavy gauge boson and SM particles [3] . 
In order to resolve the problem, the implementation of the Z 2 symmetry called T-parity to 
the model has been proposed [U E]J6] . Under the parity, almost all new particles are T-odd, 
while the SM particles are T-evemr H Thanks to the symmetry, dangerous interactions stated 
above are prohibited [7j. Furthermore, the lightest T-odd particle (LTP) becomes stable, 
which is electrically and color neutral, and has a mass of (9(100) GeV in many little Higgs 
models with T-parity @]. Therefore, these models provide a good candidate for dark matter 
[80. 

In this article, we study signatures of the littlest Higgs model with T-parity (LHT) [5j [6] 
at the Large Hadron Collider (LHC), which is expected to explore various new-physics models 
[TU| ITT] . The LHT is the simplest model realizing the LH scenario with the T-parity, and 
considered to be an attractive reference model. Since the LHC is a hadron collider, new 
colored particles have an important role to explore physics beyond the SM. As shown in the 
next section, top-partners are necessarily introduced in the LH models, which are responsible 
for the cancellation of quadratically divergent corrections to the Higgs boson mass from top 

^One important exception is the top-partner T+, which is a T-even new particle as shown in the next 
section. 

# 2 For UV completion of T-parity models, see [5]. 
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loop diagrams. Furthermore, masses of these partners are expected to be less than ~ 1 
TeV, and the partners will be copiously produced at the LHC [12]. Therefore, we consider 
the productions of the top partners at the LHC with a realistic simulation study, and show 
that these signatures are clearly distinguishable from SM backgrounds. Furthermore, we find 
that it is also possible to test the LHT by investigating a non-trivial relation among the 
signatures. We also consider how accurately model parameters of the LHT are determined, 
and discuss its implication to the property of the LTP dark matter such as how precisely the 
relic abundance of the dark matter is estimated with the LHC data. 

This paper is organized as follows. In the next section, we briefly review the littlest Higgs 
model with T-parity paying particular attention to the gauge-Higgs and top sectors of the 
model. We also present representative points used in our simulation study. Signatures of the 
LHT at the LHC are shown in Sec. [3l especially focusing on the pair production of T-even 
top partner, the single production of T-even top partner, and the pair production of T-odd 
partner. The test of the LHT is discussed in Sec. HI where we investigate a non-trivial relation 
among the signatures obtained in the previous section. We also discuss the implication of 
the result to the LTP dark matter phenomenology. Sec. [5] is devoted to summary. 

2 Model 

In this section, we briefly review the littlest Higgs model with T-parity focusing on gauge- 
Higgs and top sectors of the model. (For general reviews of little Higgs models and their 
phenomenological aspects, see P3HH].) We also present a few representative points used in 
our simulation study at the end of this section. 

2.1 Gauge-Higgs sector 

The littlest Higgs model with T-parity is based on a non-linear sigma model describing an 
SU(5)/SO(5) symmetry breaking. The non-linear sigma field E is given as 

S = e 2m// E , (2.1) 

where / ~ 0(1) TeV is the vacuum expectation value of the breaking. The Nambu-Goldstone 
(NG) boson matrix n and the direction of the breaking Eo are 

/ H/y/2 $ \ /0 1\ 

n = HyV2 H T /V2 , Eo = 1 . (2.2) 
V $ f H*/y/2 / \1 0/ 

Here, we omit the would-be NG fields in the n matrix. An [SU(2) xU(l)] 2 subgroup in the 
SU(5) global symmetry is gauged, which is broken down to the diagonal subgroup identified 
with the SM gauge group SU(2)ixU(l)y. Due to the presence of the gauge interactions 
and Yukawa interactions introduced in the next subsection, the SU(5) global symmetry is 
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not exact, and particles in the II matrix become pseudo NG bosons. Fourteen (= 24 — 10) 
NG bosons are decomposed into representations 1 © 3 © 2 ±1 / 2 © 3±i under the electroweak 
gauge group. The first two representations are real, and become longitudinal components of 
heavy gauge bosons when the [SU(2) xU(l)] 2 is broken down to the SM gauge group. The 
other scalars 2 ±1 / 2 and 3±i are a complex doublet identified with the SM Higgs field (H in 
Eq. (12.21) ) and a complex triplet Higgs field ($ in Eq. (12.21) ). respectively. 
The kinetic term of the £ field is given as 



f 2 



^S-iv / 2{^(WS + SW T )+^'(BS + SB T )} , (2.3) 



where W = W?Q* (B = BjYj) is the SU(2) j (U(l) j ) gauge field and g (g') is the SU(2) L 
(U(l)y) gauge coupling constant. With the Pauli matrix a a , the generator Qj and the hyper- 
charge Yj are given as 

, (o a 0\ 

F 1 =diag(3,3,-2,-2,-2)/10, (2.4) 




Y 2 = diag(2,2,2,-3,-3)/10. (2.5) 
It turns out that the Lagrangian in Eq. (12.31) is invariant under the T-parity, 

n <-> -nnn, w? «-> b x «-> b 2 , (2.6) 

where Q = diag(l, 1, —1, 1, 1). 

This model contains four kinds of gauge fields. The linear combinations W a = (W* + 
W 2 a )/\/2 and B = (Bi + B 2 )/V2 correspond to the SM gauge bosons for the SU(2) L and 
U(l)y symmetries. The other linear combinations Wjj = (W® — W£)/v2 and Bh = (Bi — 
B2)/V2 are additional gauge bosons, which acquire masses of O(f) through the SU(5)/SO(5) 
symmetry breaking. After the electroweak symmetry breaking with (H) = (0,v/ \^2) T , the 
neutral components of and Bh are mixed with each other and form mass eigenstates Ah 
and Z H , 



Z H \ = fcos9 H -sm6 H \ (W% 
A H ) cos6 H j \B H 



(2.7) 



The mixing angle 8h is given as 



,2 



tan** = * mi2 ^ = ~ -0.15^, (2.8) 

mn - + V ( m n _ m 22) + 4mf 2 J 

where m n = g 2 f\c) + 7)/8, m 12 = gg' j\\ - c})/8, m 22 = g' 2 f 2 (5c} + 3)/40, and c f = 
cos(v / 2t>/ /). Since the mixing angle is considerably suppressed, A H is dominantly composed 
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of Bh- Masses of gauge bosons are 





(2.11) 



(2.10) 



(2.9) 



m u + m 22 + \ (mn - «i 22 ) 2 + 4mf. 
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(2.12) 



mn + m 22 - y (m n - m 22 ) 2 + 4m 2 



12 



0.2^' 2 / 2 



(2.13) 



As expected from the definitions of A H , Z H , and Wh, the new heavy gauge bosons behave 
as T-odd particles, while SM gauge bosons are T-even. 

A potential term for H and $ fields is radiatively generated as [TJ [8] 



Main contributions to u 2 come from logarithmic divergent corrections at 1-loop level and 
quadratically divergent corrections at 2- loop level. As a result, /i 2 is expected to be smaller 
than f 2 . The triplet Higgs mass term, on the other hand, receives quadratically divergent 
corrections at 1-loop level, and therefore is proportional to f 2 . The quartic coupling A is 
determined by the 1-loop effective potential from gauge and top sectors. Since both fi and A 
depend on parameters at the cutoff scale A ~ 47r/, we treat them as free parameters in this 
paper. The mass of the triplet Higgs boson $ is given by m| = A/ 2 = 2m 2 l f 2 /v 2 , where nih 
is the mass of the SM Higgs boson. The triplet Higgs boson is T-odd, while the SM Higgs is 
T-even. 

Gauge-Higgs sector of the LHT is composed of the kinetic term of £ field in Eq. (12.30 and 
the potential term in Eq. (I2.14p in addition to appropriate kinetic terms of gauge fields W°", 
Bj and gluon G. It can be seen that the heavy photon Ah is considerably lighter than other 
T-odd particles. Since the stability of A H is guaranteed by the conservation of T-parity, it 
becomes a good candidate for dark matter. 

2.2 Top sector 

To implement T-parity, two SU(2) doublets and and one singlet ur are introduced 
for each SM fermion. Furthermore, two vector-like sing lets C/M and are also introduced 
in the top sector in order to cancel large radiative corrections to the Higgs mass term. The 
quantum numbers of the particles in the top sector under the [SU(2) x U(l)] 2 gauge symmetry 
are shown in Table [TJ All particles are triplets under the SM SU(3) C (color) symmetry. 



V{H, $) = A/ 2 Tr _ f/H^H + - (H^H) 2 + ■■■ . 



(2.14) 
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(2, 1/30; 1,2/15) 


q (2) 


(1,2/15; 2, 1/30) 




(1,8/15; 1,2/15) 




(1,2/15; 1,8/15) 




(1,8/15; 1,2/15) 


u m 

R 


(1,2/15; 1,8/15) 


UR 


(1,1/3; 1,1/3) 







Table 1: Quantum number for [SU(2) x U(l)} 2 for particles in the top sector. 



With these particles, Yukawa interactions which are invariant under gauge symmetries 
and T-parity turn out to be 



Ai/ 
2^2 



(Q (2) S ) t t ]X t ky - Qf } E ix E J u R - X 2 f U ( L ] U R n) + h.c, (2.15) 



n=l 



where = (q^,U^\o) T , 



-a \u 



L i°L ) > 



and S = Eq^S'^Eq. The indices z, j, 



run from 1 to 3, while x, y = 4, 5. The coupling constant Ai is introduced to generate the top 
Yukawa coupling and X 2 f gives the vector-like mass of the singlet U^ n \ Under T-parity, q( n ' 
and transform as q^ <-> —q^ and —U^ 2 \ thus T-parity eigenstates are given as 



q 



(±) 



(g« =F g < a >) 
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(±) 



rr(l) -prr(2) 



1 

V2 v ' 1 ' ; ' ^ L(R) ~ 7! 

In terms of the eigenstates, mass terms in Eq. (12.151) are written as 



-Ai 



fU 



u R -X 2 f (u { L +) U R +) + U 



(-)rr(-) 
L U R 



+ h.c. 



(2.16) 



(2.17) 



T-even states u + and U+ form the following mass eigenstates 



t L 



cos (3 — sin (3 
sin /3 cos (3 




t R 
T + r 



cos a — sin a 
sin a cos a 




(2.18) 



Mixing angles a, /3 and mass eigenvalues mt, rriT + are given as 

2B t C t 



tana 



tan (3 



m t 



A, 



{A 2 + B 2 _ ^ 

2A t B t 



Ai/A 
A? 



2- 



A* - (4 2 " B, 



CD 



xt + xy 

A1A2 



±y/A2 + B* + C?-A t 



V2 



V2 



Bf + Ci + A t 



v/aF+a 2 



A? + A|/, 



(2.19) 



where A t = s/AJ/v^, S t = (l + c / )A 1 //2, C t = A 2 f, and A t = ((A* + B* + C*) 2 -4A*Ci)V 2 
with sj being = sin(\/2i; //). The t quark is identified with the SM top quark, and T + is 



its T-even heavy partner. On the other hand, the T-odd fermions U l- and Ur- form a Dirac 
fermion, XL, whose mass is given by = X 2 f. The remaining T-odd quark q_ acquires 
mass by introducing an additional SO(5) multiplet transforming nonlinearly under the SU(5) 
symmetry. Therefore, the mass term of the quark does not depend on Ai and A2. In this 
paper, we assume that the g_ quark is heavy enough compared to other top partners, and 
that it is irrelevant for the direct production at the LHC experiment. (For the phenomenology 
of the q_ quark, see [15].) Finally, it is worth notifying that the T-odd partner of top quark 
(XL) does not participate in the cancellation of quadratically divergent corrections to the 
Higgs mass term. The cancellation is achieved by only loop diagrams involving t and T + 
quarks. 

2.3 Representative points 

In this paper, we focus on X± productions at the LHC. For this purpose, we need to choose 
representative points to perform a numerical simulation. In order to find attractive points, 
we consider those consistent with electroweak precision measurements and the WMAP ex- 
periment for dark matter relics^ 3 . 

We consider a x 2 -function to choose representative points; 



where 0^ s , and AC^ S are experimental result, theoretical prediction, and the error 

of the observation for observable O. We consider following eight observables; W boson 
mass (m w = 80.412±0.042 GeV), weak mixing angle (sin 2 ^ e ff pt = 0.23153±0. 00016), lep- 
tonic width of the Z boson (Ti = 83.985±0.086 MeV) [TB], fine structure constant at the 
Z pole (q-'K) = 128.950±0.048), top quark mass (m t = 172.7±2.9 GeV) [T7|, Z boson 
mass (m z = 91.1876±0.0021 GeV), Fermi constant (G F = (1.16637±0.00001)xl0" 5 GeV" 2 ) 
|18j . and relic abundance of dark matter (Odm^ 2 = 0.119±0.009) [19]. On the other hand, 
theoretical predictions of these observables depend on seven model parameters; /, A2, nih, 
a~ 1 {mz)i Gf, m z, and m t . (For the detailed expressions of the theoretical predictions, see 
[3, IE]). In order to obtain the constraint on / vs. A2 plane, we minimize the \ 2 function in 
Eq. (12.201) with respect to parameters m^, a _1 (m^), Gp, mz, and m t . In other words, we 
integrate out these parameters from the probability function P = e~ x I 2 . 

The result is shown in Fig. [H where the constraints on / and A2 at 99% confidence level 
(x 2 = 11.34) are depicted. The region A 2 < 1 is not favored due to electroweak precision 
measurements, because a large mixing angle between t and X + is predicted in this region, 
which leads to a significant contribution to the custodial symmetry breaking. The region / < 

# 3 We consider the WMAP constraint only for choosing a representative point. In fact, the model does not 
have to satisfy the constraint, because, for instance, dark matter may be composed of other particles such as 
the axion. 




(2.20) 
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Figure 1: Constraints to the littlest Higgs model with T-parity on / vs. A2 plane at 99 % confidence 
level. The degree of fine-tuning to the quadratic coupling of the Higgs field is also shown as light 
shaded regions. Cross marks 1,2, and 3 are representative points for our simulation study. 
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570 GeV, which corresponds to rriA H < rriw, is not attractive because the pair annihilation 
of Ah into gauge-boson pair is kinematically forbidden. Here, we should comment on other 
parameters integrated out from the probability function. It can be easily seen that a~ 1 (mz), 
Gp, m z , and m t are almost fixed due to the precise measurements of these observables. 
Furthermore, once (/, A2) is fixed, is also fixed by the WMAP observation, because 
the annihilation cross section of dark matter is sensitive to m^. Here and hereafter, at 
each (/, A2) point, we use values of these parameters which minimize the x 2 -function. The 
degree of fine-tuning to set the Higgs mass on the electroweak scale is also shown in the 
figure. As mentioned in the previous subsections, the quadratic coupling of the Higgs field 
fi 2 is generated radiatively. One of main contributions comes from the logarithmic divergent 
correction of a top-loop diagram, which yields [20] 

2 m k A?A| ( A 2 \ , , 

^- 3 4tATT% lo H 1 + <j' (2 - 21) 

where A ~ Airf is the cutoff scale of the model. We used the ratio F = 100 x (2m 2 h )/(^) % 
to estimate the degree of fine-tuning. It can be seen that too large / and A2 are not attractive 
from the view point of the fine-tuning. 

Representative points used in our simulation study are shown in Fig. [1] and their details 
can be found in Table [2J Masses of A H and T±, cross sections for T± pair and single T + 
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Point 1 Point 2 Point 3 



/ (GeV) 
A 2 

sin (3 

m h (GeV) 

m A[I (GeV) 

m T _ (GeV) 

m T+ (GeV)_ 

a(pp -> T_f_ + X) (pb) 

a(pp^T + f + + X) (pb) 

a ( pp ^ T+ + X) (pb) 

a(pp^f + + X) (pb) 

Br(T+ -> 

Br(T + -> Zi) 

Br(T+ -> /it) 

Br(T + - T_A H ) 



570 


600 


570 


1.0 


1.1 


1.4 


0.20 


0.16 


0.11 


145 


131 


145 


80.1 


85.4 


80.1 


570 


660 


798 


772 


840 


914 


1.26 


0.54 


0.17 


0.21 


0.13 


0.07 


0.29 


0.15 


0.05 


0.14 


0.07 


0.02 


50.8 % 


50.8 % 


53.3 % 


21.1 % 


21.8 % 


23.6 % 


15.8 % 


17.4 % 


19.1 % 


12.3 % 


10.0 % 


4.03 % 



Table 2: Representative points used in our simulation study. 



productions, and branching ratios of T + decay are also shown in each representative point. 
Note that the T_ quark decays into the stable A H and the top quark with almost 100% 
branching ratio. 



3 Signals from the LHT Events 

Now, we consider the T + and T_ production processes and their signals at the LHC. At the 
LHC, there are two types of T + production processes, pair production and single production 
processes, both of which are important. Thus, in the following, we discuss these processes 
separately. In addition, we also discuss the T_T_ pair production. 

3.1 T+T+ pair production 

First, we discuss the T + T + pair production process. Once produced, T + decays as T + — > bW + , 
tZ, hZ, and A^Tl. Branching ratios for individual decay modes depend on the underlying 
parameters. However, in most of the cases, Br(T + — > bW + ) becomes larger than 0.5, and 
many of T + decay into bW + . Thus, in the experimental situation, the analysis using the 
decay mode T + — > bW + is statistically preferred. In such a case, the t quark production 
events become irreducible background. We will propose a set of kinematical cuts suitable for 
the elimination of background. 

For the T + T + production process, the most dangerous background is the ti production 
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which has larger cross section than the T + T + production^ 4 !. Thus, we need to develop kine- 
matical cuts to suppress the tt background. We propose to use the fact that the jets in 
the signal events are likely to be very energetic because they are from the decay of heavy 
particles (i.e., T + or T+). Consequently, the signal events are expected to have large M e s, 
which is defined by the sum of transverse momenta of high px objects and missing transverse 

, (miss) 

momentum p^, : 

M cS = Y,Pt + E Pr + E & + Pt^- (3- 1 ) 

jets leptons photons 

In our study, only the jets with p T > 30 GeV are included into the high p T objects in order 
to reduce the contamination of QCD activities. We expect that the number of background 
events can be significantly reduced once we require that M e g be large enough; in the following, 
we will see that this is indeed the case. 

Once the backgrounds are reduced, the T + T + production events are reconstructed rela- 
tively easily. Here, we concentrate on the dominant decay mode T + — > bW + . Then, the signal 
events are primarily from the process pp — > T + T + , followed by T + — > bW + and T + — > bW~. 
In particular, in order to constrain the mass of T + , we use the process in which one of the W- 
boson decays hadronically while the other decays leptonically. At the parton level, the final 
state consists of two 6-jets, two quark jets from W ± , one charged lepton and one neutrino 
from W T . Thus, the signal events are characterized by 

• Several energetic jets, 

• One isolated lepton, 

• Missing px (due to the neutrino emission). 

Using the fact that, in the signal events, the missing momentum is due to the neutrino 
emission, we reconstruct two T + systems, which we call T+^-system and T_j. had -system; 
here, the Tji^-system (T^^-system) consists of high pt objects which are expected to be 
from T + or T + whose decay is followed by the leptonic (hadronic) decay of the IV-boson. 
To determine T+ ep ^- and T^^-systems, we first assume that all the missing p? is carried 
away by the neutrino. With this assumption, the neutrino momentum p u (in particular, 
the ^-component of p v ) is calculated, requiring (p t + p v ) 2 = m^. Then, we define T^ lep ^- 
system as the charged lepton, reconstructed neutrino, and one of the three leading jets, 
while T+ had ^-system is the rest of the high px objects. Since there is a two-fold ambiguity 
in reconstructing the neutrino momentum, there exist six possibilities in classifying high-j>r 
objects into 2+ p - and T^^-systems. Using the fact that T+ ep ^- and T_[ had ' ) -systems have 
the same invariant mass in the ideal case, we choose one of the six combinations with which 
|M (ic P ) — M (had)| is minimized, where M p ep ) and M (had) are invariant masses of Ti lep ' ) - 



# 4 We use the leading order calculation of the tt production cross section which is 460 pb. 
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and Tj. -systems, respectively. The distributions of the invariant masses of T| iew - and 
J , | had ) -systems are expected to provide information about the T + mass. 

In order to demonstrate how well our procedure works, we generate the events for the 
processes pp — ► T + T + and pp — > ti (as well as those for pp — ► jT + and — > jT + ) with 
£ = 100 fb _1 . The parton-level events are generated by using the MadGraph/MadEvent 
packages [2T] , which utilizes the HELAS package [22J. Then, Pythia package [23] is used for 
the hadronization processes and the detector effects are studied by using the PGS4 package 
[21]. In order to study the T + T + pair production process followed by the decay processes 
mentioned above, we require that the events should satisfy the following properties: 

1-0: Three or more jets with px > 30 GeV, and only one isolated charged lepton. 
In addition, we adopt the following kinematical cuts: 
1-1: Pt,i > 50 GeV (with p^ t i being the transverse momentum of the charged lepton), 
1-2: M eff > 1800 GeV, 
1-3: IM^) - M^had)! < 100 GeV. 

Notice that the third cut is to eliminate combinatorial backgrounds. We found that, after 
imposing these kinematical cuts, events from the jT + and jT + production processes are 
completely eliminated. Then, we calculate the distributions of M ( ha d). The results are 

shown in Fig. [2j As one can see, the distributions have distinguishable peaks at around 
M" T (had) ~ m T+ . In addition, ti backgrounds are well below the T + T + signal. Thus, from the 

distribution of M (had), we will be able to study the properties of T+. 

One important observable from the distribution of M (had) is the mass of T + ; once we 
obtain the peak of the distribution, it will provide us an important information about itit + ■ 
To see the accuracy of the determination of m^, we consider the bin Mb; n — |AMbm — 
M (had) < Mbin + k^M hin . Then, we calculate the number of events in the bin as a function 

of the center value M^n with the width AMbi n being fixed. The peak of the distribution 
is determined by which maximizes the number of events in the bin. We applied this 
procedure for AM b i n = 20 — 60 GeV (with C = 100 fb). Results for a set of signal and 
background events are shown in Table [3j With repeating the Monte-Carlo (MC) analysis 
with independent sets of signal samples, we found that the difference between the position 
of the peak and the input value of is typically 10 — 20 GeV or smaller. Thus, we expect 
a relatively accurate measurement of rriT + ■ In discussing the test of the LHT model at the 
LHC, we quote 10 and 20 GeV as the uncertainty of my + and discuss the implication of the 
measurement of rriT + - 

3.2 Single production of T + 

As well as the pair production, the single production processes pp — > jT + and jT + have 
sizable cross sections at the LHC. (Here, j denotes light quark jets.) Such processes were 
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Figure 2: Distribution of M (had) for the Points 1 — 3 (from the top to the bottom) with C 

100 fb _1 . The shaded histo grams are the background distribution, while the solid ones are 
signal + background. 
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Point 1 Point 2 Point 3 



AM bin = 30 GeV 755 GeV 834 GeV 913 GeV 

AM bin = 40 GeV 757 GeV 839 GeV 918 GeV 

AM bin = 50 GeV 741 GeV 837 GeV 910 GeV 

AM hin = 60 GeV 745 GeV 847 GeV 912 GeV 



Table 3: Peak of the M {had) distribution for AM bin = 30, 40, 50, and 60 GeV. 



discussed in [25J in the framework of the original littlest Higgs model without the T-parity, 
which pointed out that the discovery of T + may be possible by using this process. (See also 
|26j.) Here, we reconsider the single production process for the test of the LHT model. 

So far, we have discussed that the information about the mass of T + can be obtained 
by studying the T + T + pair production. Concerning the property of T+, another important 
parameter is the mixing angle j3, which determines the interaction between T + and weak 
bosons (i.e., W ± and Z). Importantly, the cross sections for the processes pp — > jT + and 
jT + are strongly dependent on (5. In particular, since these processes are dominated by 
the t-channel W /± -boson exchange diagram (with the use of b- or 6-quark in the initial-state 
protons), the cross sections are approximately proportional to sin 2 /?. Thus, if the cross 
sections of the single production processes are measured, it provides an information about 
the mixing angle f3. Although pp — > jT + and jT + have different cross section, their event 
shapes are very similar (if we neglect the charges of high px objects). In the following, we 
consider how we can measure the total cross section a pp ^jT + + o~ pp ^jf+- 

As we have already discussed, once produced, T + dominantly decays into b and W + . 
Thus, if we consider the leptonic decay of W + , there exist two energetic quarks and one 
charged lepton (as well as neutrino) at the parton level in the final state. Since the mass 
of T + is relatively large, the fe-jet is expected to be very energetic in this case. Thus, if we 
limit ourselves to the cases with the leptonic decay of W + , the single production events are 
characterized by 

• Two (or more) jets, one of which is very energetic (due to the 6-jet), 

• One isolated lepton, 

• Missing px (due to the neutrino emission). 

As we will see, the cross section of the background events are relatively large, so it is necessary 
to find a useful cut to eliminate the backgrounds as much as possible. 

One of the possible cuts is to use the invariant mass of the u bW system. In the signal 
event, the dominant source of the missing transverse momentum is the neutrino emission 
by the decay of W + . Thus, as we have discussed in the study of the T + T + pair production 
process, we can reconstruct the momentum of neutrino (and hence that of W + ). Then, we 
can calculate the invariant mass of the bW ± system. In such a study, we presume that the 
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highest pt jet is the 6-jet because, at least at the parton level, the transverse momentum of 
the 6-quark from the decay of T + is much larger than that of the extra quark. Then, since we 
expect that the mass of T + is well understood by the study of T + T + pair production process, 
as discussed in the previous subsection, we only use the events with relevant value of the 
invariant mass to improve the signal-to-background ratio. 

To estimate how well we can determine the cross section of the single production process, 
we generate the signal and background events for £ = 100 fb _1 . In [25J, it was pointed out 
that the most serious backgrounds are from ti production process as well as from the single 
production of the top-quark. Thus, in our study, we take account of these backgrounds. 

Once the event samples are generated, we require the following event shape: 

11-0: The number of isolated lepton is 1, the number of jets (with p? > 30 GeV) is 2. 

In the next step, as in the case of the T + T + pair production, we reconstruct the momen- 
tum of the neutrino assuming that the transverse momentum of the neutrino is given by the 
observed missing p?. In reconstructing the neutrino momentum p u , there exists two-fold am- 
biguity; we denote the reconstructed neutrino momenta p„ (i = 1,2). For each reconstructed 
momentum, we calculate the invariant mass of the bW system: 

Mi% = ^(p S i+Pi + ^y, (3-2) 

postulating that the highest pr jet corresponds to the b jet. Even though one of is with 
the wrong p„ , we found that, in the signal event, the typical difference between M 6 ^ and 

(2) 

are relatively small compared to that in the background events. Thus, we reject the 
events unless \M^ — is small enough. 

We also comment on another useful cut to eliminate the ti background. In the ti back- 
ground events, the highest pt jet is likely to be from the overlapping of several hadronic 
objects from different partons if the pt is required to be very large. In our analysis, the cone 
algorithm (with AR = 0.5) is used to identify isolated jets. Then, if several partons from the 
decay of top quark or IV-boson are emitted in almost the same direction, hadronized objects 
from those partons are grouped into a single jet, which may be identified as the 6-originated 
jet in the present analysis. One of the method to reject such a background is to use the 
jet-mass variable, which is the invariant mass of the jet constructed from all the (observed) 
energy and momentum that are contained in the jet. The jet-mass of such a jet is likely to 
be much larger than that of the 6-jet. As we will show, the number of background from the 
ti production process is significantly reduced if the jet mass is required to be small enough. 

Now, we show the results of our MC analysis. In our analysis, we use the following 
kinematical cuts: 

II-l: p T ,i > 100 GeV, p { T iss) > 100 GeV, 

II-2: PT,ji > 300 GeV, and Mj 1+ j 2 > 500 GeV, with Mj 1+ j 2 being the invariant mass of total 
jets, 
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II-3: Mji < 50 GeV, with Mji being the jet mass of the leading jet, 
II-4: \M^-M^\< 50 GeV. 

In Fig. [31 we plot the distribution of the "averaged" invariant mass of the bW system: 

M bW = \ (M« + Mg) . (3.3) 

As one can see, the distribution from the signal events is peaked at around Mf,w ~ ^t+, while 
the background distribution is rather flat. In addition, at around M bW ~ m T+ , the number 
of signal events becomes significantly larger than that of background in particular when the 
parameter sin (3 is relatively large. In such a case, the number of the single production events 



can be extracted from the distribution by using, for example, the side-band methocr 5 . 

In Table HI with the data for the Point 2, we show the number of events in the event 
region, which we define rriT + — 50 GeV < Mbw < + 50 GeV, and those in the sidebands, 
m T+ - 150 GeV < M bW < m T+ - 50 GeV and m T+ + 50 GeV < M bW < m T+ + 150 GeV, 
after imposing the kinematical cuts mentioned above. Assuming that the numbers of signal 
and background events in the signal region are determined by using the sideband events, and 
that the cross section for the single T + production process can be obtained from the number 
of events in the signal region, the single T + production cross section may be determined with 
the uncertainty of 10 — 20 %. (The uncertainty here is statistical only.) 

Using the result of the determination with the T + T + pair production process, the 
information about the cross section can be converted to that of the mixing angle (3. If the 
uncertainties in the theoretical calculation of the cross sections are under control, we obtain 
a constraint on (3. Since the cross section for the single production process is proportional to 



sin 2 p, sin f3 is determined with the accuracy of 5 — 10 % if the cross section is determined 



with the accuracy of 10 — 20 %F 6 l In the next section, we discuss the implication of the 
determination of (3 at this level in testing the LHT model. 

Before closing this subsection, we comment on the uncertainties which we have neglected 
so far. As we have mentioned, the single production process occurs by using the b or b in the 
sea quark of the initial-state proton. Thus, for the theoretical calculation of the cross sections, 
it is necessary to understand the parton distribution functions for the b and b quarks (as well 
as those of lighter quarks). Information about the parton distributions of the 6-quark may 
be obtained by using the single top (and anti-top) productions. As we have seen, significant 
amount of single top productions occur at the LHC (which has been seen to be one of the 
dominant backgrounds to the single T + production process). Since the single top production 
also occurs by using the b quark in proton, information about the parton distribution function 



# 5 It should be also possible to constrain the mass of T + from the peak of the distribution of M\,w ■ I n this 
paper, we will not discuss such a possibility. 

# 6 The cross section also depend on the mass of T + . Thus, the constraint on the cross section should provide 
a constraint on the (3 vs. rriT + plane. In our discussion, for simplicity, we only consider the f3 dependence 
of the cross section by using the fact that the mass of T+ is expected to be determined from the T+T+ pair 
production process. 
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Figure 3: Distribution of M&w for the Points 1 — 3 (from the top to the bottom). The shaded 
histograms are the background distribution, while the solid ones are for signal + background. 



15 



Lower Sideband Event Region Upper Sideband 











Signal 


ti 


jt + jt 


Signal 


ti 


Jt + ji 


Signal 


it 


jt + ji 


II-O 








313 


21706 


13509 


522 


12585 


8609 


116 


7810 


5362 


II-O, 


1 






108 


3366 


376 


234 


2352 


363 


44 


1747 


237 


II-O, 


1, 


2 




45 


428 


53 


144 


446 


76 


14 


440 


86 


II-O, 


1. 


2. 


3 


30 


30 


47 


114 


27 


50 


8 


21 


69 


II-O, 


1, 


2. 


3,4 


21 


12 


18 


84 


11 


12 


2 


3 


16 



Table 4: Number of the signal events /tt background/single top background events in the event 
region (my + — 50 GeV < M b w < niT + + 50 GeV) as well as in the lower and upper sidebands 
(m T+ - 150 GeV < M bW < m T+ - 50 GeV and m T+ + 50 GeV < M bW < m T+ + 150 GeV, 
respectively). Point 2, where niT+ = 840 GeV, is used. 



of b will be obtained by studying the single top production process. In this paper, we do not 
go into the detail of such study, but we just assume that the parton distribution function 
of b will become available with small uncertainty once the LHC experiment will start. We 
also note here that it is also important to understand the efficiency to accept the single 
production events (as well as the background events) after the cuts, whose uncertainties have 
been neglected in our discussion. 

3.3 T_T_ pair production 

For the study of the LHT model at the LHC, it is also relevant to consider the T-odd top 
partner, T_, and the lightest T-odd particle, Ah- For the study of T-odd particles, it is 
important to consider the T_T_ pair production process, which was discussed in [271 128] . 
Here, we reconsider the importance of this process for the test of the LHT model. 

At the LHC, T_ is pair produced via pp — > T_T_, then decays as T_ — > tA H . Since A H 
is undetectable, the T_ production events always result in missing p? events and hence the 
direct measurements of the masses of T_ and Ah are difficult. 

One powerful method to study rriT_ and rriA H is the so-called Mti analysis [29] , combined 
with the hemisphere analysis [30] . If the t and i systems are somehow reconstructed, one can 
constrain and vtia h from the distribution of the so-called Mti variable. For the event 
pp — > T_T_ followed by T_ — > tAn and T_ — > iAff, the M^ 2 variable is defined as 

M| 2 (m A J= min [max {m£( P * t , p^; m Ajl ), M^q*, , q^; m A j}l , (3.4) 

where the transverse mass Mt is defined as 



M T (p T , p£*;77U H ) = a/(| Pt | 2 + m?)(|p^| 2 + m\ H ) - p*,p£* (3.5) 

with rfiA H being the postulated mass of Ah to calculate Mr 2 . In the above expression, 
Pt = (PxiPyify an d Qt = (Qxj Qyj 0) are transverse momenta of t and i, respectively, which 
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are obtained from the reconstructed top systems. The reconstruction of the top systems is 
possible with sizable efficiency by using the hemisphere method [28j. In addition, and 
q^ H are postulated transverse momenta of the final-state Ah particles, which satisfy 

pt + q* T + pt h + <4 H = o. (3.6) 

In the calculation of Mt2, Pt h an d are varied under the above constraint to minimize 
the quantity in the square bracket of Eq. (13.41) . 

The important property of the Mt2 variable is that, if rfiA H is equal to rriA H , the upper 



end-point of the distribution of M?2 is given by tti t W\ Thus, once many samples of T_T_ 
production events become available at the LHC, it will be possible to determine the distri- 
bution of the Mt2 variable for each value of rfiA H - The distribution of the M^ 2 variable for 
the T_T_ production process was studied in [28J with the choice of rfiA H = ttiAh- In our 
discussion, we use the Mt2 analysis to constrain rriA H and rriT_, so it is necessary to study 
the distribution of the Mt2 variable for various values of m^. 

To see how the distribution of the Mt2 variable depends on rriA H , we generate the TLTL 
events (as well as ti backgrounds) and derive the distribution of Mt2- Here, we intend to use 
the events: 

pp -> T_T_ -> tA H tA H -> bW + bW~A H A H -> bqq'bq"q"'A H A H} (3.7) 
and we adopt the kinematical cuts used in [2"g] : 
III-O: No isolated leptons, 
III-l: # iss) > 200 GeV, and ^ miss) > 0.2M efr . 

Notice that large missing px is expected due to the emission of two Ah particles. Then, in 
order to reconstruct two top systems, we use the hemisphere analysis with which all the high 
Pt objects are assigned to one of two hemispheres, Hi and H 2 , so that 



d{p h 17 Pi) < d(p H 2,Pi) ■ for Vz G Hi 
d(pH 2 ,Pi) < dipHnPi) ■ for \/iEH 2 ' 

where pn I is the momentum of the J-th hemisphere which is defined as 



(3.8) 



# 7 For a general value of rfiA H , the upper end-point of the Mt2 distribution is given by 



2 

M(max) / ~ \ T- 
T2 { m A H ) = 



2mT_ \ y 2vht_ 

This can be used to check the validity of the MC analysis. 



m a — mr. 
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and 



M \ {E Hl -\PH I \cos6 Ii )E Hl 

d(PH I ,Pi) = 7^ — , (3.10) 

with On being the angle between p# 7 and p«. (For the details to construct the hemispheres, 
see [27].) In the following, the first hemisphere is defined as the one which contains the 
leading jet. Once two hemispheres are determined, we impose the following cuts to eliminate 
backgrounds: 

III-2: Numbers of jets (with px > 30 GeV) in Hi and H 2 are either equal to or smaller than 
3. 

III-3: pr,/fj > 200 GeV (J = 1, 2), where Pt,Hj is the transverse momentum of the hemisphere 
Hi. 

III-4: 50 GeV < M Hl < 190 GeV (I = 1,2), where M Hl is the invariant mass of the 7-th 
hemisphere (i.e., Mhj = \Jphj)- 

As shown in [2B], with the cuts III-O — III-3, peaks around ~ m t are obtained in the distri- 
butions of Mm and M# 2 - Then, postulating that the momenta of t and t are given by those 
of two hemispheres, we calculate the distribution of the Mti variable defined in Eq. (13. 4p 
for several values of the postulated mass ifiA H - Here, we use the underlying parameters for 
the Point 2. The results for rfiA H = 0, 100 GeV, and 200 GeV, for which the theoretically 
expected end-points are 648 GeV, 664 GeV, and 708 GeV, respectively, are shown in Fig. 
SJ Here, the distributions shown in the figure include contributions from the ti background; 
however, we have checked that there is no contamination of the ti events at the end-point 
region. 

As one can see, the position of the upper end-point changes consistently with the the- 
oretical value of the end-point. Thus, by using the My 2 variable, we expect to obtain a 
constraint on the w,a h vs. plane, which can be transferred to a constraint on the A2 vs. 
/ plane. In order to derive the constraint, it is necessary to understand how well the position 
of the upper end-point can be determined. Detailed properties of the distribution of the Mt2 
variable should depend on the kinematical cuts as well as on the detector performances. An 
extensive study of the fitting function to determine the end-point is beyond the scope of this 
paper. Here, we simply use the quadratic function to estimate the end-point. For example, 
for ffiA H = 100 GeV (for which the theoretical prediction of the end-point is 664 GeV), the 
end-point is estimated as M^ ax) = (664 ± 9) GeV (M^ ax) = (676 ± 3) GeV) using the data 
with 550 GeV < M T2 < 650 GeV (580 GeV < M T2 < 680 GeV). Thus, in the following 
discussion, we adopt the error of 10 — 20GeV in the determination of the end-point, although 
a better result may be possible if a detailed study of the shape of the end-point is performed. 



18 



- 300 



C 
> 



50 



i i i i i i i i i i i i i i i i i i i i i i i i i i i i i 




Point 2 




500 



i i i i i i i i i i i i i i 



1111 



550 600 650 700 750 800 
M T2 (GeV) 



Figure 4: Distribution of the Mt2 variable for friA H = (darkly shaded: blue), 100 GeV (lightly- 
shaded: pink), and 200 GeV (solid line). 



4 Test of the LHT Model 

Now we discuss how and how well we can test the LHT model using the results obtained 
in the previous section. As we discussed in Section [2} the LHT model is parametrized by 
two parameters, / and A2. Thus, if there exists three or more observables, a non-trivial test 
becomes possible. 

In the following, we adopt the Point 2 as the underlying parameter point, and assume 
that m,T + , sin f3, and the end-point of the Mj- 2 variable can be experimentally determined as 

[ m T+] Point 2 

±5m T+ , (4.1) 

[sin (3] p i n t 2 ±5 sin/?, (4.2) 

[M^% oult2 ±5M^\ (4.3) 

where [• ■ -Jpoint 2 denotes the value in Point 2. From the discussion in the previous section, 
we adopt the following uncertainties of the quantities mentioned above: 

• Case 1: 



m T+ = 
sin (3 = 



5m T+ = 20 GeV, (4.4) 
5 sin (3/ sin (3 = 10%, (4.5) 
SM^ x) = 20 GeV, (4.6) 
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Figure 5: Expected constraints on the / vs. A2 plane for Cases 1 and 2 (upper and lower, respec- 
tively). Point 2 is used as the underlying parameter point. Constraints from the measurements of 
rriT + , sin (3, and M^™ ax) are given by the lightly-shaded (pink) region, region between the dashed 
lines, and darkly-shaded (blue) region, respectively. The star in the figure is the underlying point. 



• Case 2: 

Sm T+ = 10 GeV, (4.7) 

S sin p/ sin/? = 5%, (4.8) 

5M^ x) = 10 GeV, (4.9) 

In the Case 2, smaller uncertainties are adopted compared to the Case 1. 

In Fig. [5J we show the allowed region on the / vs. A2 plane for the Cases 1 and 2. As 
one can see, measurements of mr + , sin/3, and provide three different constraints on 
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the / vs. A2 plane. It should be noticed that, because each of the constraints gives a narrow 
band on the / vs. A2 plane, we can quantitatively test if the observed signals are consistent 
with the predictions of the LHT model; if the three bands meet at a single point, as shown 
in Fig. El it gives a quantitative confirmation of the LHT model. 

It is also notable that the measurements of rriT + , sin P, and M^ ax) give accurate determi- 
nations of / and A2. For example, reading the lower and upper bounds on these parameters 
from the allowed region in the Case 1 (Case 2), we obtain the constraints 566 GeV < f < 
624 GeV and 1.03 < A 2 < 1.20 (584 GeV < / < 613 GeV and 1.06 < A 2 < 1.15). One of the 
implications is that, with the determination of /, we can also determine rriA H in the LHT 
model. (See Eq. ( 12. 13ft .) Since Ah is a very weakly interacting particle, the direct determi- 
nation of its mass is difficult as discussed in the previous section. Thus, the determination 
of / gives an important information about ra/^ n . 

Finally, we discuss an implication to cosmology. Ah is a viable candidate of dark matter. 
The thermal relic density of An strongly depends on the pair annihilation cross section of 
An] in the present case, An pair-annihilates into weak boson pair via the s-channel exchange 
of the Higgs boson. The pair annihilation cross section is obtained once / and are known. 
As we have discussed, / can be determined with the studies of the top partners. In addition, 
at the LHC, it is expected that the Higgs boson will be found and its properties will be 
studied in detail. For example, if = 130 — 150 GeV, the Higgs mass will be determined 



with the uncertainty of ~ 200 MeV [TQj, [TT] in the following, we assume that the Higgs 



mass can be determined with the accuracy of 200 MeV at the Point 2. Then, combining 
the information about the top-partners and the Higgs boson from the LHC, it will become 
possible to reconstruct the thermal relic density of An- Comparison of the theoretically 
calculated relic density and observed dark matter density provides an important test of the 
cosmological scenario in the framework of the LHT model; if the theoretical prediction of the 
relic density is consistent with the dark matter density observed, it will be a strong indication 
of the scenario where Ah is dark matter. 

To see how well we can perform this test, we calculate the thermal relic density O^; 
the contours of constant flA H h 2 (with h being the Hubble constant in units of km/sec/Mpc) 
are shown in Fig. [6] on / vs. plane. When / < 570 GeV, Ah becomes lighter than W^. 
In such a case, the pair annihilation cross section of Ah is extremely suppressed, resulting 
in very large value of VLA H h 2 . On the contrary, for / > 570 GeV, flA H h 2 is found to have 
mild dependence on / and m^. In the same figure, we also show the expected constraints 
on / and m^. As one can see, determination of the / parameter plays an important role in 
reconstructing the dark matter density. In particular, we can see that, combined with the 
precise measurement of the Higgs mass, QA H h 2 can be reconstructed very accurately in the 
Case 2 where the masses of top partners and mixing parameter j3 are well determined; with 
the determination of rrih and / for the Case 2 mentioned above, the density parameter is 
constrained to be 0.118 < VIa h h 2 < 0.126. (The underlying value of VIa h h 2 is 0.120.) On the 
contrary, in the Case 1 where the uncertainty in / is relatively large, bound on the density 



# 8 For a discussion of Higgs phenomenology in the LH models, see 
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Figure 6: Contours of constant Q<A H h 2 on / vs. plane. Contours are for Q,A H h 2 = 0.06, 0.08, 
0.10, 0.12, 0.14, 0.16, 0.18, and 0.20 from above. Expected bound on / is shown in the shaded 
region; lightly-shaded (pink) region is for the Case 1 while the darkly-shaded (blue) region is for 
the Case 2. The dotted lines are the expected upper and lower bounds on the Higgs mass. 



parameter is found to be flA H h 2 > 0.118. Thus, in such a case, flA H h 2 cannot be bounded 
from above. This is mainly due to the fact that we chose the underlying value of rriA H close 
to mw] with larger value of rriA H , a better reconstruction of flA H h 2 is expected even with a 
larger uncertainty in /. 

5 Summary 

In this paper, we have studied the T + T + pair, single-T + and T_T_ pair productions at the 
LHC in the framework of the littlest Higgs model with T-parity, by performing a numerical 
simulation on three representative points. For T + T + pair production process, the main SM 
background comes from ti production. We have developed kinematical cuts to suppress the 
ti background, and found that the signal events can be well extracted from the background. 
We have shown that an accurate determination of the mass of T + is possible. For single-T + 
production, we have also proposed a set of kinematical cuts to suppress the SM backgrounds 
which are from ti pair production and single-t production, and shown that the signal events 
can be well reconstructed. From the measurement of the single- T + production cross section as 
well as the measurement of rriT + in the T + T + pair production, we can obtain the information 
on the mixing parameter (sin/3) between T + and top quark. For T_T_ pair production, 
studying the upper end-point of the Mt2 distribution (M^T^), a certain relation between 
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rriA H and ttit_ is obtained. 

Since the top sector in the LHT is parametrized by two parameters, / and A2, each 
measurement of these three observables provides a relation between / and A2. We have 
shown that the measurements of the three observables give non-trivial determinations of the 
parameters / and A2, and hence a quantitative test of the LHT model can be performed at 
the LHC. 

In the LHT model, Ah is a viable dark matter candidate. Since the thermal relic density 
of Ah strongly depends on the pair annihilation cross section of Ah into weak boson pair via 
s-channel exchange of the Higgs boson, the masses of Ah and Higgs boson are important to 
calculate the thermal relic density of Ah- Using the facts that Ah mass can be determined 
by the parameter /, and that not only the discovery of Higgs boson but also the measurement 
of the Higgs mass are expected at the LHC, we have shown that the relic density of Ah can 
be calculated very accurately by using the LHC results. This will provide an important test 
of the cosmological scenario where Ah becomes dark matter. 

Our studies here suggest not only that the LHC has a great potential to discover the heavy 
partner of top quark which is responsible for the cancellation of the main quadratically 
divergent contribution to the Higgs mass parameter, but also that the LHC can provide 
important measurements of the observables that would lead us to a crucial tests of the LHT 
model. 
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